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O ■ Abstract 

m 

' In their 1978 paper, Dyson, Lieb, and Simon (DLS) proved the existence of 

Neel order at positive temperature for the spin-S" Heisenberg antiferromagnet on 
the d-dimensional hypercubic lattice when either S > 1 and d > 3 or S" = 1/2 
and d is sufficiently large. This was the first proof of spontaneous breaking of a 
continuous symmetry in a quantum model at finite temperature. Since then the 
ideas of DLS have been extended and adapted to a variety of other problems. In 
this paper I will present an overview of the most important developments in the 
study of the Heisenberg model and related quantum lattice systems since 1978, 



o 

'i> • including but not restricted to those directly related to the paper by DLS. 
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1 Introduction 

The purpose of this article is to review a number of important developments in the 
mathematical physics of quantum spin systems that have occured since Barry Simon 
first got involved with the subject. A comprehensive account of such developments 
clearly would not fit in the space alotted here. Therefore, I have selected seven topics 
that either have a relation to the seminal paper by Dyson, Lieb and Simon jT^] or 
that for some other reason I believe interesting and important enough to deserve to be 
discussed in some detail here. In Section El I list a number of works I regretfully had to 
leave out. 

There is one development that stands out by its absence: a proof of long-range order 
in the Heisenberg ferromagnet at positive temperature in dimensions > 3 has eluded 
the valiant attempts of many. One can trace of some of these attempts in the literature 
jl4| ISSl I2n] , and the authors of [T^ also comment on a erroneous proof they thought 
they had, but certainly most failed approaches have not left a written record. It is not a 
hopeless problem, in my opinion. It is much like other open phase transition problems 
such as hard-core bosons at finite density and temperature in M^, or the Xy-model in a 
uniform external field. Proving long-range order in the Heisenberg ferromagnet should 
be simpler than the problem of proving ferromagnetic long-range order in an itinerant 
electron model such as the Hubbard model [HSI iS2j- There are of course many other 
open problems in quantum spin systems. I will mention some of them in the following 
sections. And many interesting problems have been solved in the almost three decades 
since DLS1978. The ferromagnetic phase transition is a nagging reminder that we lack 
the mathematical arguments to tackle some of the most basic and clearly formulated 
problems of equilibrium statistical mechanics. Other problems in mathematical physics 
suffer from being physically not well-understood, from being mathematically ill-posed, 
from being too complex, but not the Heisenberg ferromagnet. 

In the years following DLS, Barry Simon made several other contributions to statis- 
tical mechanics that will stand the test of time: elegant estimates, useful inequalities, 
and a great book [TUl [7H1 EOl EZI EH] • Is there a chance Barry will once more turn 
his attention to statistical mechanics? The few years in the late seventies and the early 
eighties he worked in statistical mechanics were a spectacular success. Maybe I can en- 
tice him with the following spectral problem for the Heisenberg ferromagnet which lies 
somewhere in the middle between his more recent interests and full-fledged statistical 
mechanics. 
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The problem is to prove that the Hamiltonian of the spin- 1/2 isotropic Heisenberg 
ferromagnet on 'L^ in the GNS representation of one of its translation invariant pure 
ground states (which are the fully polarized states) has purely absolutely continuous 
spectrum except for the ground state eigenvalue 0. In one dimension this result follows 
from the work of Babbitt and Thomas [Hj- So the open problem is for > 2. To 
formulate the problem precisely we introduce the relevant Hilbert space as P{J-''^), 
where JF"' is the set of finite subsets of Z'^. Let {fix \ X e JF*^} be the canonical 
orthonormal basis of Ti. One should think of Vlx as the state vector with all spins 
pointing up except for those finitely many that are located at x G X, which are down. 
The Hamiltonian H is then the unique selfadjoint operator with domain containing the 
canonical orthormal basis and such that 

Hn^ = 0, Hflx = - ^xvj{y}\{.}) , for all ^ X G . 

x&X ,y^X ,\x—y\=l 

Note that, up to an overall sign, H is also the generator of the symmetric simple 
exclusion process on Z''. Another expression for H is given in ()2.1|) . It is easy to 
see that H is non-negative definite and from the definition it is obvious that is an 
eigenvalue. It is not hard to prove that the eigenvalue is simple. The conjecture is 
that for all d > 1 the positive real line is absolutely continuous spectrum of H . In 
particular, there are no strictly positive eigenvalues. 

2 General properties 

The mathematical setting for studying quantum spin systems rigorously, as described 
in textbooks such as [75j; and (79j, was developed in the late sixties starting with 
Robinson [71j. This is not to say that there are no rigorous results about quantum 
spin models before that date. In particular Lieb and coworkers wrote several landmark 
papers that continue to be important today (see, e.g., [171110]). In this general formalism 
it was possible to obtain several fundamental theorems for large classes of quantum 
spin models. An early example is the uniqueness of the equilibrium state at any finite 
temperature for one-dimensional models with short-range interactions by Araki [7j. In 
this section we first briefly review the main elements of the general setup and then 
discuss some recent results in this context. 

For the purpose of our discussion, a quantum spin is any quantum system with a 
finite-dimensional Hilbert space of states. A quantum spin system consists of a finite 
or infinite number of spins, which we will label by the elements of a set V . When 
V is an infinite set, typically corresponding to the vertices of a lattice or a graph, 
one often considers families of quantum spin systems, labeled by the finite subsets 
X d V . Certain properties are easily stated for infinite sets V directly, but to define 
the dynamics and specific models we will first assume that V is finite. In this case, the 
Hilbert space of states is 

= (g)C"^ 

x&V 
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where the dimensions Ux > 2 are related to the magnitude of the spins by Hx = 2sx + 1, 
and Sx G {1/2, 1, 3/2,...}. For each spin, the basic observables are the complex n x n 
matrices, which we will denote by Mn- The algebra of observables for the system is 
then 

xev 

Given a Hamiltonian, a self-adjoint observable Hy = Hy G Ay^ the Heisenberg dynam- 
ics of the system is defined as follows: for any t G M, an automorphism at is defined on 
A by the formula 

a'dA) = u;r{t)AUv{t) 

where Uv{t) = e~**^^. 

One of the most important examples of a quantum spin system is the Heisenberg 
model. In general, this model is defined on a graph (K, E) which consists of a set of 
vertices V and a set of edges E. E is a, set of pairs of vertices denoted by e = {xy) for 
x,y & V. Let 5*, i = 1,2,3, denote the standard spin Sx matrices associated with the 
vertex x, and for each edge e = {xy), let Jxy G M be a coupling constant corresponding 
to e. The Heisenberg Hamiltonian (also called the XXX Hamiltonian) is then given by 

Hy = — ^ JxySx ■ Sy , (2.1) 
ixy)eE 

where Sx denotes the vector with components Sl,S^,S^. The most commonly studied 
models are those defined on a lattice, such as Z'^, with translation invariance. For such 
Hamiltonians, the magnitude of the spins is usually constant, i.e., Sx = s, the edges 
are the pairs (xy) such that |a; — |/| = 1, and J^y = J. Depending of the sign of J, the 
Heisenberg model is said to be ferromagnetic (J > 0) or the antiferromagnetic (J < 0). 

For extended systems, i.e., those corresponding to sets V of infinite cardinahty, more 
care is needed in defining the quantities mentioned above. Hamiltonians are introduced 
as a sum of local terms described by an interaction, a map $ from the set of finite 
subsets of V to Ay, with the property that for each finite X G V, $(X) G Ax and 
^{X) — $(X)*. Given an interaction the Hamiltonian is defined by 

Hy=Y^ HX). 
XCV 

Infinite systems arc most often analyzed by considering families of finite systems, in- 
dexed by the subsets of V, and taking the appropriate limits. For example, the C*- 
algebra of observables. A, is defined to be the norm completion of the union of the local 
observable algebras [jxcv^^- 

To express the decay of interactions and correlation functions we need a distance 
function on V. Let V be equipped with a metric d. For typical examples, V will be a 
graph and d will be chosen as the graph distance: d{x, y) is the length of the shortest 
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path (least number of edges) connecting x and y. The diameter, D{X), of a finite 
subset X C V is 

D{X) = max{d{x, y) \ x,y E X}. 

In order for the finite-volume dynamics to converge to a strongly continuous one- 
parameter group of automorphisms on .4, one needs to impose a decay condition on 
the interaction. For the sake of brevity, we will merely introduce the norm on the 
interactions that will later appear in the statement of the results described in detail 
here. For weaker conditions which ensure existence of the dynamics see fT^ 17^1 17^ EH] • 
We will assume that the dimensions are bounded: 

N = supn^^ < oo, 

x&V 

and that there exists a A > such that the following quantity is finite: 

< oo. 

For translation invariant systems defined on a lattice with an interaction $ such that 
II $11 A < oo for some A > 0, Lieb and Robinson 53^ proved a quasi-locality property of 
the dynamics, in the sense that, up to exponentially small corrections, there is a finite 
speed of propagation. More precisely, the obtained an estimate for commutators of the 
form 

where t G M, A G Ax, B G Ay, and X,Y dV. Clearly, such commutators vanish 
if t = and X fl F = 0. Quasi-locality, or finite group-velocity, as the property is 
also called, means that the commutator remains small up to a time proportional to the 
distance between X and Y . We will now formulate an extension of the Lieb-Robinson 
result to systems without translation invariance or even an underlying lattice structure. 
It will be useful to consider the following quantity 

\\[at{A),B]\\ 

CB{x,t) ■= sup r— , 

A£A^ ll^ll 

for xdVjtE'R, Be Ay- The basic result is the following theorem |64j . 
Theorem 2.1. For x G V", t G M, and B G Ay, we have the bound 

CB{x,t) < e2|*lll*II^CB(x,0) 

+ J2 e-^"(^'^)(e2|*lll*IU-l)CB(y,0). 

Our proof avoids the use of the Fourier transform which seemed essential in the 
work by Lieb and Robinson and appeared to be the main obstacle to generalize the 
result to non-lattice (y,d). 
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If the supports of A and B overlap, then the trivial bound || [tj (A), 5] || < 2||y4|| ||i?|| 
is better. Observe that for B G Ay, one has that Cb(2/, 0) < 2||i?|| Xriv), where xy 
is the characteristic function of Y, and therefore ii x ^ Y, then one obtains for any 
A G Ax a bound of the form 

\\[Tt{A),B]\\ < 2\Y\ \\A\\\\B\\ (e^ 1*1 11*11^ - l) g-^'^^^'^), 

where 

d{x, Y) = mm{d{x, y) \ y E Y} . 
Moreover, for general local observables A G Ax, one may estmate 

\\[n{A),B]\\<N'\''\\\A\\J2CB{x,t), 

in which case. Theorem 12.11 provides a related bound. 

Theorem 12 . II has been further extended to provide bounds in the case of long-range 
interactions by Hastings and Koma pfB|. 

The next general result relies on the Lieb- Robinson bound to prove that a nonvanish- 
ing spectral gap above the ground state implies exponential decay of spatial correlations 
in the ground state. This result can be regarded as a non-relativistic analogue of the 
Exponential Clustering Theorem in relativistic quantum field theory j21j. The idea that 
a Lieb- Robinson bound can be used as a replacement for strict locality in the relativistic 
context is natural but, as far as I know, it was used in the literature for the first time 
only recently by Hastings in [S3] . 

In the physics literature the term massive ground state implies two properties: a 
spectral gap above the ground state energy and exponential decay of spatial correlations. 
It has long been believed that the first implies the second, and the next theorem proves 
that this is indeed the case. The converse, that exponential decay must be necessarily 
accompanied by a gap is not true in general. Exceptions to the latter have been known 
for some time and it is not hard to imagine that a spectral gap can close without 
affecting the ground state. 

For simplicity of the presentation, we will restrict ourselves to the case where we 
have a representation of the system (e.g, the GNS representation) in which the model 
has a unique ground state. This includes most cases with a spontaneously broken 
discrete symmetry. Specifically, we will assume that our system is represented on a 
Hilbert space 7i, with a corresponding Hamiltonian H > 0, and that f2 G 7i is, up to a 
phase, the unique normalized vector state for which HQ = 0. We say that the system 
has a spectral gap if there exists 6 > such that spec(if) fl (0, 6) = 0, and in this case, 
the spectral gap, 7, is defined by 

7 = sup{(5 > I spec(if) n (0, 5) = 0}. 

Our theorem on exponential clustering derives a bound for ground state correlations 
which take the form 

{n,Anf,{B)n) (2.2) 
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where b >0 and A and B are local observables. The case 6 = is the standard (equal- 
time) correlation function. It is convenient to also assume a minimum site spacing 
among the vertices: 

inf d{x,y) =: a > 0. (2.3) 



We proved the following theorem in 

Theorem 2.2 (Exponential Clustering). There exists fi > such that for any 
X ^ y & V and all A e Ax, B G Ay for which BVt) = 0, and b sufficiently small, 
there is a constant c{A, B) such that 



\{n,Aai,{B)^)\<c{A,B)e 

One can choose 



'Ud(x,y)( 1+ T^^L^T^ 



4||$||a + 7 

and the hound is valid for < 76 < 2^d{x, y). 

The constant c(A, 5), which can also be made explicit, depends only on the norms 
of A and 5, (in its more general form) the size of their supports, and the system's 
minimum vertex spacing a. For 6 = 0, Theorem 12.21 mav be restated as 

I (fi, AB^) - (fi, Afi) (fi, 5fi) I < c(A, B)e-^'^^'''y\ (2.6) 

Note that there is a trivial bound for large 6 > 

\{Q, Aaib{B)Q)\ < \\A\\ \\B\\e-^^. (2.7) 

In the small 6 > regime, the estimate (|2.4p can be viewed as a perturbation of (12. 
Often, the important observation is that the decay estimate fl2.4|) is uniform in the 
imaginary time ib, for b in some interval whose length, however, depends on d{x,y). 

In a recent work, Hastings and Koma have obtained an analogous result for models 
with long range interactions [36 . 

It has been known for some time that the spontaneous breaking of a continuous 
symmetry, such as the SU{2) rotation symmetry of the XXX Heisenberg model ()2.H) . 
precludes the existence of a spectral gap in the infinite volume ground state. This 
is often called the Goldstone Theorem of statistical mechanics |45| IHHj . Matsui 
adapted the proof of the Goldstone theorem to show the absence of a spectral gap 
above the Ill-interface ground states of the XXZ ferromagnet defined in (j3.8|) in all 
dimensions d > 2, a. result anticipated by Koma and Nachtergaele who had previously 
proved it for = 2 

In dimensions < 2 spontaneous continuous symmetry breaking cannot occur at 
positive temperatures. This is called the Mermin- Wagner- Hohenberg theorem a general 
version of which was proved by Frohlich and Pfister [^H] and for systems of fractal 
dimension c? < 2 by Koma and Tasaki In the next section we review some results 
on existence of spontaneous symmetry breaking in quantum spin models at T = and 
T > 0. 
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3 Neel order and other applications of reflection 
positivity 

One of the main results in the 1978 paper by Dyson, Lieb, and Simon (DLS), and 
an extension by Jordao Neves and Fernando Perez [20] to include ground states, and 
refinements by Kennedy, Lieb, and Shastry [1^11101; is the existence of long range order 
in the isotropic spin-S" antiferromagnet on hypercubic lattice if" at sufficiently small 
positive temperatures in dimensions > 3, and any spin magnitude > 1/2, and also 
in the ground state (T = 0), if d > 2 and S" > 1, or d = 3 and S = 1/2. 

DLS was a breakthrough based on another breakthrough, namely, the method of 
infrared bounds to prove continuous symmetry breaking developed for classical systems 
by Frohlich, Simon, and Spencer [211 ISOl ■ The method was soon generalized to cover a 
wide variety of models [2H1 12ZI that have a property called refiection positivity. There is 
no a priori reason why the infrared bound philosophy should require refiection positivity, 
a property not shared by the Heisenberg ferromagnet jSTj. But in spite of serious 
attempts to find ways to prove infrared bounds for quantum spin systems not relying 
on refiection positivity, no such method has been found to date, and we still do not 
have a proof of long range order in the ferromagnetic Heisenberg model at positive 
temperature. 

For the XXZ model with Hamiltonian 



an Ising-type phase transition has been proved by Frohlich and Lieb [23; for d > 2 
and sufficiently large A (since their proof assumed refiection positivity). Kennedy 
P7] showed how to prove the phase transition for all A > 1 without using refiection 
positivity but by instead relying on a version of the Peierls argument using a clever 
cluster expansion method. Thus it may appear that to make further progress one should 
circumvent the refiection positivity property. This is certainly true in part. It is worth 
mentioning however, that the use of refiection positivity in quantum lattice models 
pioneered in DLS has been successfully adapted to solve a number of other problems. 
Let us just mention a few interesting examples: dimerization in the Hubbard- Peierls 
model on the ring and the spin-Peierls problem [22], the fiux-phase problem [^ 130] , 
and singlet ground states of quantum lattice systems [IHl 1221 1211 122] 

4 The classical limit and beyond 

In a very elegant paper, [IH], Lieb proved the existence of the classical limit for a 
broad class of quantum spin systems for which the Hamiltonians can be considered as 
a multilinear function, H, of the standard spin operators. The main result of Lieb's is 
then upper and lower bounds for the quantum partition function with the spin operators 




(3.8) 
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taken to be of magnitude 5*. To state the result we define 

where the spin operators in the Hamiltonian are taken to be of magnitude S and the 
trace of the corresponding Hilbert space of dimension (25* + 1)'^', and where V is fixed. 
The corresponding classical partition function is defined by 

ZciP) = I e-'^^^^i" n 

where dQx is the normalized invariant measure on the 2-sphere (unit vectors in M^). 
Then, Lieb proved 

ZciP) < Zq{P, S) < Zcm + |)) (4.9) 

with c a universal constant. This result suffices to prove convergence of the mean 
free energy (assuming it exists for the classical system), as well as convergence of the 
expectation of averaged quantities (intensive observabels) . 

In jTB] Simon generalizes Lieb's result by proving the same relations for systems 
of G-spins, where G is an arbitrary compact Lie group. The Hamiltonian is assumed 
to be a multilinear function operators X*, which form the basis for the Lie-algebra 
of G in a fundamental representation of weight A, copied at each x & V. He then 
considers the sequence of quantum systems where each X* is replaced by 7r(X*), with 
TT a representation of weight nX, ra = 1, 2, 3, . . .. Simon's main result is then 

ZciP) < ZQiP,n) < + ^)) (4.10) 

where Zc is again a classical partition function for a suitable classical system which he 
defines. 

The bounds fl4.9j) and ()4.1()|1 by themselves are not good enough to deduce a phase 
transition for the quantum model at sufficiently large 5* or n, assuming that the classical 
model has a phase transition, although one may expect such a result. Very recently, 
Biskup, Chayes, and Starr P] refined the Berezin-Lieb inequalities, which are at the 
core of the results of Lieb and Simon, to include estimates of matrix elements e~^^ in 
coherent state vectors, not just for the partition function. This enables them to prove 
the existence of phase transitions in the quantum model for sufficiently large spin, 
assuming that both the classical and the quantum model are refiection positive and 
that chessboard estimates prove a phase transition for the the classical model. Another 
refinement of the classical limit, very much akin to studying the central limit theorem 
as a refinement of the law of large numbers, was undertaken in j6n[ l6T] . 

5 Spin chains and Haldane's conjecture 

Quantum spin chains, i.e., systems defined on the one- dimensional lattice Z, and their 
ground states are a fascinating subject by themselves. They have received a lot of 
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attention in recent years due to their relevance for quantum information theory. But 
even before quantum information theory became again fashionable, the subtleties and 
great variety in behavior attracted the attention of condensed matter physicists and 
mathematical physicists alike. Haldane gave the subject a big impetus when he pre- 
dicted that the Heisenberg antiferromagnetic chain has qualitatively different ground 
states depending on whether the magnitude of the spins is integer or not The 
Bethe Asatz solution of the spin- 1/2 chain shows a unique ground state with polyno- 
mial decay of correlations and no spectral gap. Haldane conjectured that the spin-1 
chain would also have a unique ground state (in the thermodynamic limit), but with a 
spectral gap above it and exponentially decaying correlations (this is sometimes called 
the Haldane phase). Affleck and Lieb [HI proved that for half-integer spin chains a 
unique ground state indeed imphes the existence of low-lying excitations. A related 
implication was proved by Aizenman en Nachtergaele |3|: a unique ground state must 
necessarily have no faster than polynomial decay of correlations. 

That the Haldane phase exists was rigorously demonstrated by Affleck, Kennedy, 
Lieb, and Tasaki, who proposed the spin-1 chain with Hamiltonian (now called the 
AKLT model) 



and proved that it has a unique ground state, which they explicitly constructed, with 
exponential decay and a spectral gap above the ground state 

Inspired by the AKLT model and a construction by Accardi PP, Fannes, Nachter- 
gaele, and Werner introduced Finitely Correlated States [21] and proved that they 
provide the exact ground state of a large family of spin chains, including the AKLT 
model, with similar properties. 

All correlations of Finitely Correlated States (also called Quantum Markov States or 
Matrix Product States) can be explicitly computed by a transfer matrix type formula. 
Let Az be the algebra of local observables of the spin chain with each spin represented 
by a copy of ^ = M„, and let B be another finite-dimensional C*-algebra, say B = 
Mk- Let K : A ® B ^ B, he a completely positive map that is unit preserving i.e., 
IE(]ly4 (8) lie) = Ig. It can be shown that for such E there exists at least one state p on 
B such that 



Suppose p is such a state. It is useful to introduce the maps : B — » i3, E^(i?) = 
E(y4 ® -B). It is then straightforward to verify that, with these objects and properties, 
the following formula defines a translation invariant state lj on Az'- 



We say that u is the finitely correlated state generated by the triple {B, E, p) . 



A particulary important class of maps E are the so-called pure maps. Let V : 

® be linear and isometric, i.e., V*V = 1. Then, E{A ® B) = V* A ® BV, 




X 



(5.11) 



p(E{1a ®B))= p{B), for a\\B eB 



Lu{Ai ® ^2 O • • ■ O A^) = Tr pE^, o E^^ o ■ ■ ■ o Ea^W ■ 
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defines a pure completely positive unital map in the sense that there is no non-trivial 
decomposition E = Ei + E2, with Ei and E2 completely positive maps. Such finitely 
correlated states are called purely generated. Pure finitely correlated states can be 
generated by a pure map E 

We conclude this section by specifying the triple {B, E, p) that generates the ground 
state of the AKLT model with Hamiltonian (j5.11|) . Since it is a spin-1 chain, A = M3. 
It turns out that B = M2 and E(X) = V*XV, with V : ^ given by 

V\l/2) = ^173 |1) ® 1-1/2) - v/273|0) ® |l/2) 
r 1-1/2) = 1-1) ® IV2) - 727310) ® 1-1/2) 

where |m) denotes an eigenbasis of 5"^. V is also the, up to a phase factor, unique 
isometry satisfying the following intertwining relation 

y^(l/2)^ ^(1)^^(1/2)^ 

where D^^^ is the spin-S" representation of SU{2). The state p is given by 

6 Perturbation Theory 

A major goal in perturbation theory of quantum spin systems is to show that the 
set of interactions for which the model has a unique ground state with a non-vanishing 
spectral gap above it (in the thermodynamic limit), is open in a suitable topology on the 
space of interactions. Significant steps toward this goal have been made by a number 
of authors and the results they obtained have a scope and generality that justifies 
calling them quantum Pirogov-Sinai theory El QHj • It should be noted, however, 
that the applicability of these results is restricted to models that can be regarded as 
quantum perturbations of classical models. In some examples suitable perturbation 
methods have also been used to show the absence of gaps or eigenvalues [T7| EHl- The 
remarkable paper by Kennedy and Tasaki |3Tj was the first to make a serious attempt 
to get away from perturbing classical models. A new result by Yarotsky, which we 
discuss below, can be seen as taking that line of approach one step further. 

Yarotsky's result makes it possible to prove stability of the massive phase provided 
that there is a nearby Generalized Valence Bond Solid model |21 EI 112] that can be 
used as a reference point for the perturbation in the space of interactions. In particular, 
Yarotsky [SZ] proves that the AKLT chain of 1)5.111) is contained in an open set of 
interactions with this property. 

The general theorem proved by Yarotsky can be stated for the following class of 
models defined on Z*^, d > 1. For these models the Hilbert space T~ix 5 ^t x ^ Z , 
is allowed to be infinite-dimensional. Let Hy = 'S>xev^^^ ^'^^ finite V C Z'^, 
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be the Hilbert space associated with V. The unperturbed model has finite-volume 
Hamiltonians of the form 

H^= Yl (6.12) 

where hx is a selfadjoint operator acting non-trivially only on Hvo+x, for some finite Vq. 
The main assumption is then that there exists ^ ^ T^x such that Vl^ = 
is the unique zero-energy ground state of Hy, with a spectral gap of magnitude at least 
I Vol above the ground state. Explicitly: 

H^>0, i/W = 0, Hv>\Vo\{l-\n''y){n''y\). (6.13) 
The perturbed Hamiltonians are assumed to be of the form 

Hy = H'y+ ^^x^ + ^'x\ (6.14) 

x,Vo+xCV 

where (px^ and (px^ are selfadjoint operators on Hvo+x satisfying 

\{i^,m)\<a\\hl/'^\\\ ||0f||</?, (6.15) 

for all ijj G Dom(/iy^), and suitable constants a and (3. One can call 0*-^'' a "purely 
relatively bounded" perturbation, while (f)^''^ is simply a bounded perturbation. 

Theorem 6.1 (Yarotsky [,87j ) . Let Hy of the form \6.14^ , satisfying assumption 
\6.1^) . For all K > 1 there exists 5 = S{K,d,Vo) > such that if condition i6.1^) is 
satisfied with some a G (0, 1), and j3 = 5(1 — 0;)'''^'^+^^ then 

1) Hy has a non- degenerate gapped ground state VLy '■ Hy^v = Ey^v, o,'>T'd for some 
7 > 0, independent ofV, we have 

Hy > Ey\ny){Qy\ + {Ey + -f) {1 ~ \Qy) {QyD . (6.16) 

2) There exists a thermodynamic weak*-limit of the ground states Qy : 

{Any,ny) ^^^uj{A), (6.17) 
where A is a bounded local observable. 

3) There is an exponential decay of correlations in the infinite volume ground state uo: 
for some positive c and /i, and Ai G B{HyJ, 

\uj{AiA2) - io{Ai)uj{A2)\ <cl^^l+l^^le-'^''(^i'^^)pi||p2||. (6.18) 

4 ) If, within the allowed range of perturbations, the terms (px ( or the resolvents {hx+(j)x~ 
z)'^ in the case of unbounded perturbations) depend analytically on some parameters, 
then the ground state uo is also weak* analytic in these parameters (i.e. for any local 
observable A its expectation uj{A) is analytic). 
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Application of this result to the AKLT model (jS.llll yields the following theorem. 



Theorem 6.2. Let $ = $* g ^[o,r]- Then there exists Aq > 0, such that for all 
A, |A| < Aq, the spin chain with Hamiltonian 



has a unique infinite-volume ground state with a spectral gap and exponential decay of 
correlations. Here ^ ^[x,x+r], is $ translated by x. 

To prove this theorem, Yarotsky shows that the AKLT model itself can be regarded 
as a perturbation of a particular model, one he explicitly constructs, to which Theorem 
16.11 can be applied. 

7 Estimates of the spectral gap 

The most essential result in the AKLT paper is the proof of a non- vanishing spectral 
gap in the thermodynamic limit. Also in [21] one of the core results is an explicit lower 
bound for the spectral gap of the Valence Bond Solid models, i.e., the generalizations 
of the AKLT model introduced in that paper. A further generalization to models with 
discrete symmetry breaking and hence a finite number of distinct infinite- volume ground 
states, is studied in In that paper lower bounds for the spectral gap are obtained 
by an adaptation of the so-called martingale method of Lu and Yau jSH] to the quantum 
context. In the following we will prove a simple version of this method for obtaining 
lower bounds for the spectral gap in extended systems. Then, we will discuss some 
generalizations and improvements. 

Consider a quantum spin model on with local Hamiltonians of the form 



where h^^x+i is the translation of a fixed nearest neighbor interaction /io,i, acting non- 
trivially only at the pair of sites {x, x + 1}, with /lo.i a non-negative definite element of 
^[0,1]- Suppose that w is a zero-energy ground state of this model in the sense that 



This is the situation encountered for VBS models and there are other interesting ex- 
amples [021 uni • 

Our aim is to show that, under some rather general conditions on u;, this model has 
a non- vanishing spectral gap above the ground state uo, i.e., there exists 7 > 0, such 
that 




X 



b-l 




x=a 



(^(^x.x+i) = 0, for all X E Z. 




(7.19) 
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where is the GNS Hamiltonian of the model in the state u, which satisfies H^Vt^ = 0. 

The martingale method provides lower bounds for the smallest non-zero eigenvalue 
of the finite- volume Hamiltonians H[a,b], uniform in the volume. We start with a lemma 
that implies that this is enough to establish ()7.19p . 

Lemma 7.1. Define 

7 = sup{(5 > I spec(i7^) n (0, S) = 0}, 
assuming that the set on the RHS is non-empty and put 7 = otherwise. Then 

7 > liminf Ai(r2), (7.20) 

n>2 

where Xi{n) is the smallest non-zero eigenvalue of H[i^n]- 

Proof. Let 70 be defined by the RHS of ()7.20|) . Then, we need to show that 

for all X G ^loc- We may assume that X G A[a^b], for some a < b E Z. By using the 
fact that Hi^7c{X)Qi^ ± kerif^^, we see that 

{n^, Tx{X*)Hl7,{X)n^) = hm hm lim (fi^, 7r(X*)7r {[H^.,, [Ha„ [Hj,,,X]]]) n^) 

A3 — >£a A2 — >^ Ai — >/j 

= {n^,n{X*)Hf^_,^,^,^niX)n^). 

The last expression is the expectation of an element of ^[a-3,b+3] in uj. Let p[a~3,b+3] be 
the density matrix of the restriction of u to A[a-3,b+3] ■ Then, 

= Tr (-f/'[a-3,fe+3]P[a-3,b+3]-^*-f^[a-3,b+3]) -f^[a-3,6+3] 
> Ai(6 - a + 5) Tt p[a-3,b+3]X*H[a-3,b+3]^ 

= \i{b-a + 5){n^,7i{X*)H'^7i{X)Q^) 

As X G for all A that contain [a, b], we can take the liminf of the last inequality to 
conclude the proof. □ 

With this lemma, the problem of finding a lower bound for the spectral gap of the 
GNS Hamiltonian is reduced to estimating the smallest non-zero eigenvalue of the local 
Hamiltonians. Due to the translation invariance it is sufficient to consider intervals of 
the form [l,i^]. As before, 

L-l 



H[i,L] = ''^hx,x+i (7-21) 



x=l 
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where h^^x+i is a translation of hi^2, acting non-trivially only at the nearest neighbour 
pair {x, X + 1}, and we assume that hi^2 > and that ker Hii^] 7^ {0}. We will denote 
by 72 the smallest nonzero eigenvalue of hi^2, i-e., the gap of -f^[i,2]- It is obvious that 



L-l 



kerif[i,L] = Q kei hx,x+i (7.22) 



x=l 



For an arbitrary subset A let Ga be the orthogonal projection onto 

ker ^ hx,x+i (7.23) 

x,{x,x+l}cA 

For intervals [a,b], 1 < a < b < L, G[a^b] is the orthogonal projection onto the zero eigen- 
vectors of Yl^x=^a^^,^+i^ ^{x} ~ ^ ^- Fi'OKi these definitions it immediately 
follows that the orthogonal projections Ga satisfy the following properties: 

Ga.Ga, = Ga.Ga, = Ga, if Ai C A2 (7.24-a) 
GaiGa^ = Ga^Gai if Ai n A2 = (7.24-b) 
hx,x+i > 72(]1 - ^[^,^+1]) (7.24-c) 

Define operators En, 1 < n < L, on by 

( 1 - G[i,2] if n = 1 

En = I G[i,„] - G[i,„+i] if 2 < n < L - 1 (7.25) 
[G[i,l] if n = L 

One can then easily verify, using the properties ()7.24-a|) - ()7.24-cp . that {En \ 1 < n < L} 
is a family of mutually orthogonal projections summing up to 11, i.e.: 

L 

En = En, EnEm = Sni,nEn, ^ -^n = 1 (7.26) 

n=l 

Next, we make a non-trivial assumption. 

Assumption 7.2. There exists a constant e, < e < l/\/2, such that for all 1 < n < 
L-l 

EnG 

[n,n+l] En < e^En (7.27) 

or, equivalently, 

||G[„^„,+i]i?„|| < e (7.28) 

Note that, due to (j7.24-aj) . G[„^„+i]£'„ = G[n,n+i]G[i^n] — G[i^n+i]- This relates As- 
sumption [7i21 with Lemma 6.2 in j5T], where an estimate for ||G[n,„+i]G[i,„] — G[i,„+i]|| 
is given for general Valence Bond Solid chains with a unique infinite volume ground 
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state. The same observation also implies that [G[n,n+i]iG[i^n^ = [G[n,n+i], En], which, if 
()7.28p holds, is bounded above in norm by 2e. 

The next theorem is a special case of Theorem 2.1 in ^j. Just like Theorem 6.4 
in [21] it provides a lower bound on the gap of the finite volume Hamiltonians, but it 
achieves this in a slightly more efficient way. We will repeat the proof for the particular 
case stated here, because it is simple, short, and instructive. 

Theorem 7.3. With the definitions of above and under Assumption ]?.^ the following 
estimate holds for the smallest non-zero eigenvalue of i^[i,L] 

Ai(L) > 72(1 - V2ef (7.29) 

i.e., the spectrum of H[i Lj has a gap of at /easi 72(1 — v^e)^ above the lowest eigenvalue, 
which is 0. 

Proof. From the properties ()7.26j) of the En and the assumption that G[i^l]4' = 0, it 
immediately follows that 

L-l 

= Y,\\Enijf (7.30) 

n=l 

One can estimate in terms of {ip \ hn^n+ii^) as follows. First insert G[n,n+i] and 

the resolution {Em}: 

L-l 

IlKV'ir = (V- I (1 - Gln,n+l])Enij) + | (7.31) 

m=l 

Using (|7.24-a,j) and (!7.24-b;) one easily veryfies that E^ commutes with Gn,n+i if either 
m<n — 2 or m>n + l. In these cases EmG[n,n+i]En = Gin,n+i]EmEn = 0, because the 
En form an orthogonal family. By this observation we obtain the following estimate. 
For any choice of constants Ci, C2 > 0: 

||E„^f = (V- I (1 - G[n,n+l])En^) + ((^n-1 + K)^ | G[„,„+l] ^nV') 

< -L(^|(]i_G'[„,„+i])^) + |(^|E„^) (7.32) 

+ I EnGln,n+l]Eni^) + ^(^ | (^n-1 + ^n)'^) 

/C2 ^ 

where we have applied the inequality 

K0l|02)|<^||0lf + |||02f , 

for any c > 0, to both terms of 1)7.311) . The first term in the right side of inequality 
()7.32|) can be estimated with the interaction using ()7.24-c|) . The third term can be 
estimated with ()7.27|) . It then follows that 

(2 - Ci - -)\\En^f - C2UEn-l + S„)^f < | hn,n+l^) 

C2 Ci72 
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The term containing En-i is absent for n = 1, and = if = L. We now sum 
over n and use ()7.30|) to obtain 

(2 - ci - ^ - 2c2)\m' < -^{^ I H[,,L]^) 

C2 Ci72 

Finally put ci = 1 — e\/2 and C2 = e/ and one obtains the estimate ()7.29|) stated in 
the theorem. □ 



The theorem above, which is stated in the form given in joz] , can be generalized in 
different directions. First, one can formulate the conditions for local Hamiltonians on a 
increasing sequence of finite volumes that are not necessarily intervals, or that grow to 
infinity in a more restricted way. It is also not necessary that the model is translation 
invariant, as long as the inequality of Assumption 17. 21 holds and 72 is replaced with the 
minimum of the gaps of H^n^n+i], over all n = 1, . . . , L — 1. An application where this 
is relevant was given in |15j . 

A nice improvement of the lower bounds provided by Theorem 17. HI was made by 
Spitzer and Starr in IH21. We will state their main result without proof. It is useful to 
generalize slightly, as done in |H2], by introducing a family em,n, for m,n > 1. To this 
end, we consider finite systems on intervals of the form [a, b], with a < —m, b > n, and 
define e(m, n), by 

e{m,n) = sup{(V^, G[o,„]^) | ^ G G[^rn,o]^ = tp, G[^.m,n]i' = 0, Wi'W = 1}. 

Then, e^.n, for m, > 1, is defined by 

€m,n = sup e(m',n). 

m'>m 

Theorem 7.4 ([82]). With the same notations as in Theorem ]? . !^ we have the bound 

>^i{L) > Ai(n + m) (^1 - 2^e„,„(l - em,n)^ ■ 

Even in the case m = n = 1, this theorem gives a strictly better lower bound than 
Theorem 17.31 To see this, note that ei^i = e^, Ai(2) = 72 and, for all < e < 1/a/2, one 
has 

l-2^/e2(l-e2) 
1 < < 2. 

(1 - V2ey - 

So, depending on the value of e, the bound of Theorem 17.41 is always better than the 
one provided by Theorem 17. 31 bv up to a factor of 2. 

By somewhat different but not unrelated methods excellent estimates were obtained 
by Caputo and Martinelli for the higher spin XXZ chains in jT3] • They proved that the 
spectral gap is bounded above and below by quantities of the form constant times S, 
supporting the conjecture that the limit lims^oo7('S')/S' exists and is strictly positive 
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8 Ordering of Energy Levels 



One easily checks that the Heisenberg Hamiltonian Hy-, defined in ()2.1|) . commutes 
with both the total spin matrices and the Casimir operator given by 

Sl^ = Y^Si,i = 1, 2, 3, and C = Sy ■ Sy 

x&V 

The eigenvalues of C are S'(S' + 1) where the parameter S G {^min, ^min + 1, . . . , ^max}, 
with S'max = s^- S IS Called the total spin, and it labels the irreducible represen- 
tations of SU{2). Let Ti^^'^ be the eigenspace corresponding to those vectors of total 
spin 5*. One can show that Ti^^"^ is an invariant subspace for the Hamiltonian Hy and 
therefore, the number 

E{Hy,S) := minspec 

is well-defined. 

Supported by partial results and some numerical calculations, we made the following 
conjecture in fo7\ . 

Conjecture 8.1 ([67]). All ferromagnetic Heisenberg models have the Ferromagnetic 
Ordering of Energy Levels (FOEL) property, meaning 

E{Hy,S) <E{Hy,S'), ifS'<S. 

In Lieb and Mattis proved ordering of energy levels for a class of Heisenberg 
models on bipartite graphs, which includes the standard antiferromagnetic Heisenberg 
model. The FOEL property mentioned above can be considered as the ferromagnetic 
counterpart. To compare, a bipartite graph G = (V, E) is a graph such that its set 
of vertices V has a partition V = AU B where A (1 B = ^ and any edge (xy) G E 
satisfies either x E A and y G -B, or x G -B and y E A. For such a graph, one considers 
Hamiltonians of the form 

H = —Hy + Ha + Hb, 

where Hy, Ha, and Hb are ferromagnetic Heisenberg Hamiltonians on the graph G, 
and arbirary graphs A and B, respectively. Let Sx = XIxgx X = A,B, V. The 

Lieb-Mattis Theorem jlUl [32] then states that 

(i) the ground state energy of H is E{H, \ Sa — Sb\) 

(ii) if |5a - < 5 < S', then E{H, S) < E{H, S'). 
One can see this property illustrated in Figure ^ 

We first obtained a proof of FOEL for the spin 1/2 chain in [HZj . In that paper we also 
prove the same result for the ferromagnetic XXZ chain with SUq{2) symmetry. Later, 
in ini] , we generalized the result to chains with arbitray values of the spin magnitudes 
Sx and coupling constants Jx,x+i > 0. In short, we have the following theorem. 

Theorem 8.2. FOEL holds for all ferromagnetic chains. 
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Figure 1: The spectrum of a ferromagnetic Heisenberg chain consisting of 5 spin-1 
spins with constant couphngs. On the horizontal axis we have plottted the eigenvalue 
of the third component of the total spin. The spectrum is off-set so that the ground 
state energy vanishes. The arrows on the right, with label S, indicate the multiplets 
of eigenvalues E{H,S), i.e., the smallest eigenvalue in the subspace of total spin S. 
The monotone ordering of the spin labels is the FOEL property. On the left, we 
have indicated the largest eigenvalues for each value of the total spin. The monotone 
ordering of their labels in the range 1, . . . , 5, is the content of the Lieb-Mattis theorem 
applied to this system. 
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The main tool in the proof is a special basis of SU(2) highest weight vectors intro- 
duced by Temperley-Lieb [84^ in the spin 1/2 case and by Frenkel and Khovanov [21] in 
the case of arbitrary spin. Some generalizations beyond the standard Heisenberg model 
have been announced [Hni IHH] • 

The FOEL property has a number of interesting consequences. The first immediate 
implication of FOEL is that the ground state energy of H is E{H, S'max), corresponding 
to the well-know fact that the ground state space coincides with the subspace of maximal 
total spin. Since there is only one multiplet of total spin S'max, FOEL also implies that 
the gap above the ground state is E{H, Sma.x — ^) — E{H, S'max)- In the case of translation 
invariant models this is the physically expected property that the lowest excitations are 
simple spin-waves. 

Another application of the FOEL property arises from the unitary equivalence of the 
Heisenberg Hamiltonian and the generator of the Symmetric Simple Exclusion Process 
(SSEP). To be precise, let G = (y,E) be any finite graph, and define fi„ to be the 
configuration space of n particles, for n = 0, 1, . . . , \V\, consisting of r] : V {0, 1}, 
with ^j.gy ^7(3;) = n. For any (xy) G E, let r^y > 0. The SSEP is the continuous time 
Markov process on fi„ which exchanges the states (whether there is a particle or not) at 
X and y with rate r^y, independently for each edge (xy). The case n = 1 is the random 
walk on G with the given rates. 

Alternatively, this process is defined by its generator on /^(i7„): 



where 77^^ is the configuration rj with the values at x and y interchanged. One verifies 
L > 0, LI = 0, and therefore L generates a Markov semigroup {e'^^}t>o, such that 



where /io is the initial probability distribution on the particle configurations. It is easy 
to show that for each n there is a unique stationary measure given by the uniform 
distribution on Q^- The relaxation time, which determines the exponential rate of 
convergence to the stationary state, is given by 1/A(n), where A(n) > is the spectral 
gap (smallest eigenvalue > 0) of L as an operator on /^(f2„). 

Aldous, based on discussions with Diaconis |6jj, made the following remarkable con- 
jecture concerning A(n): 

Conjecture 8.3 (Aldous). \{n) = A(l), for all I < n < \V\ - 1. 

Assuming the conjecture, one may determine the gap by solving the one-particle 
problem. 

To make the connection with FOEL, observe 



{xy)eE 




Qfi^n) = (c')®!^! = n 



V 
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This is the Hilbert space of a spin 1/2 model on V. An exphcit isomorphism is given 
by 

/ ^ ^ = ^fiv) \v) , 

where \ri) G H is the tensor product basis vector defined by 

Si \^) = - 1/2) \^) 

Under this isomorphism the generator, L, of the SSEP becomes the XXX Hamiltonian 
H . To see this it suffices to calculate the action of L: 

Lfiv) ^ J2(^f){v)\v) 

(xy) V 
(xy) V 

= Hi) 

where t^y is the unitary operator that interchanges the states at x and and the last 
step uses 

1/2 25*2; ■ Sy =1 txy, 

and Jxy ^"^xy- 

The number of particles, n, is a conserved quantity for SSEP. Since, S^^^ = —\V\/2 + 
n, the corresponding conserved quantity for the Heisenberg model is the third compo- 
nent of the total spin. Under the isomorphism the invariant subspace of all func- 
tions / supported on n-particle configurations is identified with the set of vectors with 
5*^ = S^ax ~ "-5 which we will denote by 7i„. The unique invariant measure of SSEP 
for n particles on V, which is the uniform measure, corresponds to the ferromagnetic 
ground state with magnetization n — \ V\/2. The eigenvalue A(n) is the spectral gap of 
Hvlun- It is then easy to see that the FOEL property implies that A(n) = A(l). Since 
we proved FOEL for chains, we also provided a new proof of Aldous' Conjecture in that 
case jSl]. 



9 Other topics 

To conclude I would like to give a sample of the many topics in quantum spin systems 
that I have omitted in the preceding sections. First, I have not discussed developments 
in exact solutions of quantum spin chains by the Bethe ansatz and related techniques, 
which is a subject all by itself. I also did not elaborate on the many recent results on 
the XXZ model related to interface and droplet states. For a short review on this topic 
see jnSl- There are also some interesting connections between quantum spin chains and 
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problems in combinatorics related to counting alternating sign matrices and packed 
loop configurations on a finite square lattice with matchings specified at the boundary 
(721 ESI EH]- Another topic of recent activity is the behavior of the dynamics and 
non-equilibrium properties of quantum spin chains [SHI ESI HOI El ■ 

Quantum spin systems are an inexhaustible source of mathematical problems, many 
of which are of physical interest. Somebody should write a book about them. 
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